Signature of gravitational radiation in afterglow of short Gamma-Ray Bursts? 
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Short Gamma- Ray Bursts (GRBs) are brief intense emission of 7— rays characterized by a duration 
shorter than 2 seconds and some, if not aU of them, are hkely powered by the coalescence of 
binary neutron stars. Sometimes the remnants may be supramassive highly-magnetized neutron 
stars (magnetars) with a typical rotational period Po ~ 1 ms and the rapid rotation prevents 
the gravitational collapse until a good fraction of the rotational energy has been lost. Plausible 
observational evidence for the short-living magnetar central engine reported in the literature is 
the X-ray plateaus followed by an abrupt decline identified in some short GRB afterglow. We 
however find that with a moderate/high radiation efficiency favored by the very dim optical afterglow 
emission, the observed duration of some X-ray plateaus are significantly shorter than that expected 
in the dipole radiation model (i.e., the magnetar collapsed considerably earlier than the prediction 
of the pure dipole radiation). We suggest that the gravitational radiation of the magnetar can 
enhance the collapse considerably and thus account for the "shortened" duration of the plateau. The 
corresponding gravitational radiation signal may be detectable for the proposed Einstein Telescope. 

PACS numbers: 98.70.Rz, 04.30-w 



Thanks to the successful performance of the Swift 
satellite, our understanding of short GRBs, a kind of 

S-ray outbursts with a duration less than two seconds 
, has been revolutionized 0. At least for some short 
GRBs, the binary-neutron-star merger model Q has been 
supported by host galaxy observations and by the non- 
detection of accompanying bright supernovae [3, [H| • The 
coalescence of double neutron stars inevitably produces 
an energetic burst of gravitational radiation, which is ex- 
pected to be one of the most promising targets for current 
and the proposed gravitational wave detectors [S] . More- 
over gravitational wave detection, in principle, is able to 
pin down the nature of the remnant of the merger, either 
a stellar mass black hole or a neutron star. Though the 
prospects are promising, so far direct detection of gravi- 
tational wave radiation is unavailable yet and people are 
looking forward to the performance of detectors like ad- 
vanced LIGO. Current knowledge of the binary neutron 
star mergers is mainly based on numerical simulations as 
well as observations of short GRBs. 

In the relativistic simulations of double neutron star 
merger with a total gravitational mass ~ 2.6 — 2.8 M© 
(typical for the ten binary systems identified in the 
Galaxy), a heavy rapidly- rotating neutron star with a 
spin period ^ \ ms is formed at least temporarily 
0, Q . Since Pq is crucial for this work, let us evaluate it 
analytically following |9l|. For simplicity the double neu- 
tron stars are assumed to have a similar mass M^^ ^ 
1.35 Mq, the orbital angular momentum of the binary 
system can be estimated as Jbinary — ^/GM^^r\j2, where 



Tt ~ 6rg is the tidal radius and r^,^^ — 2GM^^/c'^ flo| . 
The angular velocity of the newly formed neutron star 
with a moment of inertia / can be estimated by f2o ~ 
X Jhi-aaxy / 1 , assuming that a fraction (x) of the orbital 
angular momentum of the neutron star binary goes into 
the remnant keeping almost all of the total mass. The 
gravitational radiation changes the energy of the binary 
system at a rate dE/dt w 1.9 x 10^^(rt/6rg,fjg)~^ erg/s 
and the merger takes a time At ^ 1.4 ms (j't/Srg.^g)^ 
[lo| . The corresponding change of the angular momen- 
tum reads AJbinary {dE /dt)M/ J2GM^^/rl and we 
havex = 1— A Jbinary/ ^binary ^ 0.9 Q. Therefore One has 
f7o ~ 2 X 10* (M^„/1.35M^)2(//1045-3 g cm^)"!. The 



,/1.35M0)2(//lo45-3 gcm2)- 
above crude but straightforward estimate does confirm 
that Pq ^ 1 nis. For stiff equation of states favored by 
the discovery of a neutron star with mass 1 .97 ± 0.04 Mq 
[m , the rapid rotation is likely able to prevent the grav- 
itational collapse [l2| and a millisecond magnetar-like 
object [l^ may survive until a good fraction of its rota- 
tional energy has been lost. The internal energy dissi- 
pation of the Foynting-flux dominated outflow launched 
by the newly-formed millisecond magnetar before its col- 
lapse may power high energy transient and then account 
for the peculiar X-ray emission (i.e., unexpected in the 
external forward shock model) following short bursts [l^l , 
in paticular the X-ray plateau lasting ^ 100 s followed by 
an abrupt drop [l5| . In a very recent systematic analysis 
of the X-ray afterglow emission of short GRBs, Rowlin- 
son et al. [l6| found magnetar signature in about half 
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of events. These authors then estimated the physical pa- 
rameters, including the initial rotation periods Pq and the 
strength of dipole magnetic field B± with the usual as- 
sumptions that the duration (luminosity in a wide energy 
range) of the plateau is the spin-down timescale (wind lu- 
minosity) of the magnetar, i.e., 

L^,^, = 10^^ erg s-i rj{-^r{-§^n^)-\ (1) 
^ 10°cm lO^^'G 1ms 
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where Rs (Bg) is the radius (surface polar magnetic field 
strength) of the star, B± = SgSina and a is the an- 
gle between the rotational and dipole axes z is the 
redshift of the GRB, and rj is the efficiency of convert- 
ing the magnetar wind energy into late prompt radiation 
and has been taken to be ~ 1 in [l^ [3l ■ The inferred 
B± arc in the range lO^'^ — 10"'^'' Gauss, in agreement 
with the magnetar hypothesis. The inferred Pq for most 
events however are 3> 1 ms, at odds with the double 
neutron star merger model. Below for illustration we 
discuss GRB 080905A and GRB 090515A, two nearest 
short events with magnetar central engine signature in 
the X-ray lightcurves. With the X-ray plateau of GRB 
090515A at a redshift z = 0.405 0, Po is inferred to be 
- 10 ms [l|. For 080905 A at a redshift z = 0.122 [3, 
Pq is estimated to be ~ 5 ms [lB|- Evidently, such spin 
periods are considerably longer than that expected in the 
merger scenario. Moreover, with a Pq ^ 1 ms, the rota- 
tion is too slow to prevent the collapse of the magnetar, 
rending the model not self-consistent. 

The current situation can be re-expressed more 
straightforwardly as that the total energy released in 
the X-ray plateau phase is much smaller than that 
expected in the merger-formed magnetar model (~ 
27t'^I/Pq ~ 10^^ erg). Such a puzzle might be solved 
in two ways. One is that rj is as low as ~ 0.01 
(see Wu et al. 2013 in preparation), then the in- 
trinsic Pq can be 1 ms since L^^^^t^i^^/{1 + z) « 
4 X 10^^ erg 77(//10^^'^g cm^)(Po/lms)~^, which yields 
Po oc ?7~^/^ for giving L^^^^ and t^^^^/{l + z). Such a smaU 
r] however may be less likely since usually the magnetic 
energy dissipation process is expected to produce (late) 
prompt emission at a moderate /high efficiency ~ 0.5 "l^, 
that is why in the literature (e.g., [lEIi^) people adopted 
ry 1. Moreover in the low 77 case, the non-dissipated 
magnetar wind will inject into the blast wave of the GRB 
outflow, drive energetic forward shock and give rise to 
bright afterglow emission i^l- GRB 090515A likely lies 
outside of its host galaxy j5[ and one might be able to at- 
tribute the very dim optical emission to the extremely low 
circum-burst medium. However for GRB 080905A such 
an argument does not apply any longer since this burst 



lies offset from the host galaxy centre by ^ 18.5 kpc but 
still in the northern spiral arm which exhibits an older 
stellar population than the southern arm 15|. Hence 
we expect a normal circum-burst medium, with a kinetic 
energy > 10^^ erg and a very low redshift z — 0.122 the 
forward shock emission should be extremely bright. The 
observations are not in support of such a speculation. Ac- 
tually att^ 8.5 hours after the trigger, the R-band emis- 
sion was at the magnitude of 24.04 ± 0.47. At late times 
it faded continually and was not detected any longer at 
i > 36 hours after the trigger [l^. Shifting the burst 
to z = 1, at the observer's time ~ 15.2 hours after the 
trigger the R-band emission should be ~ 29"^ magnitude. 
However, as shown in Fig. 15 of [2l|, the bursts with a typ- 
ical kinetic energy 10^^ erg (most are long events) are 
far brighter. We therefore conclude that for some events 
the small i] argument is disfavored and then turn to the 
other possibility that most of the rotational energy was 
carried away via the strong gravitational wave radiation. 
Rotating magnetar with a triaxial shape has a time vary- 
ing quadrupole moment and hence radiates gravitational 
waves at a frequency / = 2 /P [l3, 22] . A magnetar loses 
rotational energy through magnetic dipole torques and 
gravitational wave radiation emission |17| 



~dE,ot/dt = Tr^R^^Blf/6c^ + 32Tr^Gll,e^ f /5c^ , 

where e — 2(/xx — Iyy)/ilxx + ^yy) is the ellipticity in 
terms of the principal moments of inertia (i.e., /). We 
would like to point out that the gravitational wave radia- 
tion of magnetars formed in coUapsars driving long GRBs 
has been extensively investigated [l^, [l^] but these works 
focus on the delectability by the future detectors and no 
signature of gravitational radiation in afterglow has been 
suggested possibly because the spin periods of such a kind 
of magnetars are hard to reliably estimate. 
For 
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the rotational energy loses mainly through the dipole 
radiation and the spin down timescale is estimated by 
eq.(I2]) otherwise the rotational energy loses mainly via 
the gravitational wave radiation and the spin down 
timescale is given by 
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For t > (1 + z)min{Tdip, t^wI; ^ considerable fraction 
of rotation energy has lost and the decelerated magne- 
tar likely collapses. The duration of the plateau can 
be estimated as tpiat ~ (1 + -z) iiiin{rdip, t^^} and 
the luminosity is given by eq.([T]). Hence the puzzle 
arising in interpreting the X-ray plateau data of GRB 
080905 A and GRB 090915 A within the merger-formed 
supramassive magnetar model can be solved as long as 



3 



Eplat 



(1 + z)t^„, requiring 



0.0095 ( 



Po 



1045.3g cm2' 



^l/2(_!_0_^2t 



Iplat 



(l + z)100s' 



'1/2 



(3) 

which is considerably larger than the maximum elastic 
quadrupole deformation of conventional neutron stars, 
but comparable to the upper limit derived for crystalline 
colour-superconducting quark matter [23]. Interestingly 
the required sizable ellipticity may be attributed to the 
deformation caused by the strong magnetic field of the 
magnetar (see [2^ 1 for a summary and for other possibili- 
ties) . If the internal magnetic field is purely poloidal and 
matches the dipolar field in the exterior, the ellipticity 
may be as large as - 0.01 for Bj^ - 5 x lO^^ G 
As shown in eq.dl]), such a large B± is too large in the 
moderate/high 77 case but possible in the low 77 0.01) 
case. Alternatively a number of studies have suggested 
that the internal magnetic field Bt of the magnetar could 
be dominated by a toroidal component with a strength 
up to 10^'' G [26|, which could induce a prolate deforma- 
tion with the ellipticity e - 0.016(St/10^'^ G)^, matching 
that found in eq.®. 

In our scenario, the gravitational wave radiation has 
a frequency / « 2000 Hz and the energy spectrum is 
dEg^/df ~ TT^Izzf, hence the characteristic gravitational 
wave amplitude can be estimated as i23|] 
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which is below (but comparable to) the advanced-LIGO 
noise curve /irms = [fShif)]^^'^ t2Iji, (where S'h(/) is the 
power spectral density of the detector noise) at such a 
high frequency, where D is the distance of the burst to 
the earth. The proposed Einstein Telescope with an ex- 
pected ft-rms « 5 X 10"^^ at / = 2000 Hz [ll] may be able 
to detect the gravitational radiation signal discussed in 
this work. Interestingly, with the Einstein Telescope the 
remanent (either a black hole or a supramassive neutron 
star) formed in double neutron star merger at a distance 
D < 1 Gpc may be directly pinned down [29l| and the 
collapse process might be identified too [3(I|. Though 
the magnetar GW signal discussed in this work may still 
be undetectable if Z) > 100 Mpc, one can however test 
our model indirectly by monitoring the electromagnetic 
counterpart to estimate the total energy injected into 
the medium since one of our main prediction is that 
such an energy should not be (considerably) higher than 
that emitted in the X-ray plateau phase. 
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